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ABSTRACT 

We present and analyze a new M' detection of the young exoplanet Pictoris b from 2008 
VLT/NaCo data at a separation of ~ 4 AU and a high signal-to-noise rereduction of L' data taken in 
December 2009. Based on our orbital analysis, the planet's orbit is viewed almost perfectly edge-on 
(i ~ 89 degrees) and has a Saturn-like semimajor axis of 9.50 AXJ^'j 3 ^. Intriguingly, the planet's 
orbit is aligned with the major axis of the outer disk (Q, ~ 31 degrees) but probably misaligned with 
the warp/inclined disk at 80 AU often cited as a signpost for the planet's existence. Our results mo- 
tivate new studies to clarify how Pic b sculpts debris disk structures and whether a second planet 
is required to explain the warp/inclined disk. 

Subject headings: stars: early-type, planetary systems, stars: individual Pictoris 
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1. INTRODUCTION 

Two decades of studies have argued that the nearby, 
12 Myr-old A-type star f3 P ictoris likely harbors a 
young planetary system (e.g.. ISmith and Terrild 119841 : 
Kala s and JewittJ 119951 : iMouillet et al.l I1997T) ; recently, 
Lagrange et al.l (|2009al [2010) detected a w 9 ± 3 Mj 
planet around this star (0 Pic b). Imaged at projected 
separations of ~ 6 AU and ~ 8 AU (November 2003 and 
October-Decemb er 2009, resp ectivel y), Pic b - along 
with HR 8799e (|Marois et al.l l201Ct ICurrie et al.l 120111 1 
- may also provide a more direct comparison to the so- 
lar system's gas giants than other directly imaged planets 
which ar e at wider separations (e.g. Foma lhaut b and HR 
8799bcd lKaias et alJl200l iMarois et alJl2008| ). Current 
studies have yet to detect the planet at projected sep- 
arati ons < 0.3" (|Lagrange et alJl2009ri iFitzgerald et al.1 
2009). Data at these smaller separations could provide 
crucial constraints on the planet's orbit. 

Imaging planets at small, < 0.3" separations re- 
quires significantly reducing quasi-static speckle noise 
and wavefront errors induced by imperfect AO cor- 
rections. Advanced observing/image processing tech- 
niques like angular differential imaging (ADI) coupled 
with PSF subtraction from a locall y optimized combi- 
nation of images (LO CI) algorithm ([Marois et al.l 12005 
ILafreniere et al.l I2007D significantly attenuate speckles 
and increase sensitivity. New i nstrumentation, such a s 
the Gemini Planet Imager (GPI Ma cintosh et al.| [2008). 
will achieve far superior wavefront control. Generally, 
these efforts focus on planet imaging in the near-IR. How- 
ever, mid-IR imaging naturally overcomes some of these 
challenges as the achievable Strehl ratio is better and the 
planet-to-star contrast is most favorable. High-contrast 
imaging with Strehl ratios > 0.9 can yield at least some 
planet detections close to the telescope diffraction limit 
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(e.g. 2 A/D for HR 8799d: [Serabvn et al.l [2010h . Be- 
cause M band imaging often achieves these high Strehl 
ratios (jMinowa et al.l 120101: Iffinz et al.l [20061 ). it may be 
a promising route for detecting very young and self lu- 
minous planets at small A/D separations, despite a much 
higher sky background. 

In this Letter, we report a detection of f3 Pic b at a 
separation of ~ 0.21" extracted from archival M' band 
VLT/NaCo data taken in November 2008. We also 
present a high signal-to-noise L 1 detectio n of Pic b from 
Decem ber 2009 data first published in Lagrang e et al.l 
(2010) . We combine thes e data with recent data from 
iBonnefov et al.l (|201lD and lQuanz et al.l (|2010fl to better 
constrain the orbit and atmosphere of (3 Pic b. 

2. OBSERVATIONS AND DATA 

Our study originates from the need to test 
the ADI/LOCI re duction pipeline first presen t ed in 
ICurrie et all (|2010D and updated in ICurrie et all (|2011D 
at separations smaller than those where the pipeline had 
previously extracted planet signals (r < 0.375"). Because 
Pic b's reported project ed separation in 2009 was « 
0.3" (Lagrange "et~atlf 2010). we chose the now publicly 
available ILagrange et al.l L! band data from December 
29, 20 09 to test our code performance. Lagrang e et al.l 
(2010) discusses the details of the L' band observations. 
The total field rotation in units of the image FWHM was 
~ 3 A/D, sufficient for using our reduction pipeline. 

Figure [T] (top- left panel) shows our processed L 1 band 
image using the LOCI algorithm in annular regions of 
250 x FWHM (Na = 250) with reference images selected 
from frames with at least 0.5 x FWHM field rotation (5 — 
0.5). The planet is easily detected and is well separated 
from residual speckle noise. The planet signal-to-noise, 
determined from the dispersion in pixel intensity values 
in concentric ann uli, is SNR ~ 21 : abou t a factor of 4-5 
greater than from ILagrange et al.l (|2010f ). 

Motivated by this success, we searched for additional 
Pic data in the VLT/NaCo archive taken in ADI mode 
between 2003 and 2009, finding a set taken on November 
11, 2008 with the L27 camera. Most of these data were 
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taken in sparse aperture masking mode in K s , L', and 
M' band s over a span of ~ 4 hour s: these data were men- 
tioned in lLagrange et al.l (|2009bD as not providing good 
constraints on the companion. However, we found ~ 13 
minutes of the M' data taken in ADI mode without aper- 
ture masking at various times in between the masking 
data. Over the course of the entire observing sequence, 
the parallactic angle changed by ~ 100 degrees, or ~ 2.4- 
3 A/D at 0.2"-0.25": sufficient for image processing with 
our pipeline. 

Basic image pr ocessing of the M' band data followed 
steps outlined in iCurrie et al.l (|2011| ) . After registering 
each image and subtracting off the smooth seeing halo, 
we Fourier filtered the data to remove residual low spa- 
tial frequency noise and masked any hitherto unidentified 
bad pixels previously lost in the seeing halo. We explored 
a range of LOCI parameter space, varying 5, Na, and the 
ratio of the radial to azimuthal lengths of the subtrac- 
tion annulus (g). Because Pic b is very lumino us in 
the mid-IR (e.g. AL' 7.7; lLagrange et al.ll20Tol and 
Section 3 of this work), we focused on "aggressive" LOCI 
settings of 5 = 0.25-0.5, N A = 200-300, and g = 0.3-1, 
which better remove residual speckle noise. 

Figure Q] (top-right panel) shows our best reduced M' 
image. Pic b is clearly detected in the southwest quad- 
rant w 0.2-0.25" from' the star (SNR ~ 6). Manually 
inspecting each image between the radial profile subtrac- 
tion and final image combination steps and examining a 
signal-to-noise map of the median-combined image also 
shows that the peak does not result from latent image 
artifacts. Slightly different settings for S, Na, and g also 
yield significant detections (bottom panels). 

3. ANALYSIS 

Our new M' band detection and high signal-to-noise 
L' band detection allow new constraints on the planet's 
orbit. To derive precise astrometry needed to investi- 
gate the planet's orbit, we adopt the N aCo plate scale 
and o rientation for the L27 camera from iBergfors et al.l 
1)20111 ): 27.1 mas/pixel and a north position angle of -0.6 
degrees. These values a re nearly identic a l to tho se for the 
L27 camera quoted bv lLagrange et al.1 (l2009aft for 2003 
NaCo data and for the S27 camera from lEhrenreich et al.l 
(2010) calibrated from Trapezium data acquired closest 
in time to the Pic data: our astrometric results does 
not leverage on which calibration we use. 

To fine tune our measurements, we correct for the pho- 
tometric and astrometric bias induced by LOCI process- 
ing by comparing the imputed fluxes and positions of 
fake point sources added to registered images with com- 
puted fluxes an d centroid positions obtained after LOCI 
processing (e.g. lLafreniere et al.l l2007t iThalmann et al.1 
120091 ICurrie et al.l 120 111 ). While we lack unsaturated 
data from this run to directly confirm the PSF shape, 
unsaturated M' data taken in prior runs such as that 
for HD 158882 (March 2007) show that the AO-assisted 
NaCo M' PSF core is axisymmetric and well reproduced 
by a simple Gaussian intensity distribution. For the V 
band data, the astrometric bias is minimal, whereas 
Pic b's measured radial separation in M' band is biased 
by about +0.5 pixels (0.013"). The position angle offsets 
for both data were minimal. 

We determine the M 1 band position to be at a sep- 
aration of r = 0.210 ± 0.027" and position angle of 



211.49 ± 1.9 degrees. The L' band position is at 0.326 
± 0.013" and 210.64 ± 1.2 degrees (Table 1). Here we 
conservatively assume an uncertainty in radial separa- 
tion of one pixel for M' band and 0.5 pixels for the 
(higher signal-to-noise) L' band data. The position an- 
gle uncertainty - determined from the dispersion in val- 
ues using different centroiding estimates (e.g. cntrd.pro 
vs. gcntrd.pro) is 0.25 pixels (0.7 masxr), or 1.2 and 
1.9 degrees for V and M\ compa r able to unce rtain- 
ties for Pic b bv lLagrange et all (|2009aL I2010T ) and 
iBonnefoy et al.l (1201 lh. A ssuming that Pic is 19.3 pc 
distant (jCrifo et al.lll997l ). the planet was at a projected 
separation of 4.05 ± 0.50 AU on November 11, 2008 and 
6.29 ± 0.25 AU on December 29, 2009. 

To determine the range of allowable orbits for 
Pic b, we follow the method described in Uanson et al.1 
(|2011[ ) used to model the orbit of the low-mass brown 
dwarf companion GJ 758 B (jThalmann et al.1 120091 : 
ICurrie et al.l l2010f). somewhat similar to earlier analy- 
ses fo r Pic b ((Lagrange et al.l 120093 iFitzgerald et al.l 
2009) . In this approach, we perform a Monte Carlo sim- 
ulation comparing the astrometry to predictions from 
randomly selected orbits, where we allow all orbital pa- 
rameters to vary. The minimum \ 2 value in our simu- 
lation is x 2 ~ 1-23. Given our data's weak constraints 
on the orbital acceleration and the degeneracies due to 
the unknown line-of-sight components of planet position 
and velocity, no single 'best' orbital solution emerges. 
Rather, the best-fitting solutions describe an extended, 
well- defined family of solutions that all match the data 
equally well. We choose a cut of \ 2 < 2.23 (\ 2 < X 2 min + 
1) to represent the family of best-fitting orbits. We also 
consider the results for a cut of x 2 < 8 (% 2 for a l-cr devi- 
ation in each of the data's degrees of freedom): a family 
of 'average- fitting' solutions. From the set of models 
satisifying this criterion, we determine the median value 
of each parameter, weighted by the ratio of the mean to 
current orbital velocity for the corresponding orbit, and 
identify the weighted 68% confidence interval about the 
median. We include astrometry from the highest signal- 
to-noise data separated in time by more than ~ 3 months 
(Table 1). 

Figures [2] and [3] displays our Monte Carlo simulation 
results. For a \ 2 < 2.23 cutoff (Figure [2]), the range of 
best-fit orbital parameters (weighted median, [weighted 
68% confidence interval]) include a p = 10.99 [8.18, 15.88] 
AU, i = 89.47 [89.19, 89.69] degrees, e = 0.12 [0.03,0.31], 
and O = 30.89 [30.57, 31.17] degrees. For nearly circular 
orbits (e < 0.1), the range in semimajor axes is much 
narrower (a p ~ 8-12 AU). If we relax our fitting criteria 
to accept models with x 2 < 8 (Figure [3]), we find a p — 
9.50 [7.80, 13.43 ] AU and e = 0.10 [0.02,0.23]. More 
importantly, the inclination and longitude of ascending 
node are still nearly single valued: i = 88.93 [88.06, 89.40] 
and n = 31.32 [30.56, 32.12]. Thus, Pic b's orbit must 
be viewed almost perfectly edge on, consistent with that 
for Pic's debris disks, with a north p osition angle of ~ 
30.8 degrees for the outer debris disk (|Kalas and Jewittl 
119951 see also Boccaletti et al. 2009) but inconsistent 
with the inner disk position angle, which is offset b y ~ 5 
degrees (|Heap et al.l l2000: Gol imowski et al.ll2006h . 

Unfortunately, the M' band observations of Pic were 
taken with the star saturated within ~ 3 pixels (~ 0.6 
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FWHM) of the centroid and there were likewise no un- 
saturated standard star observations. We derive a very 
crude magnitude estimate by scaling the M' PSF of HD 
158882 to the unsaturated portion of the f3 Pic PSF and 
use HD 158882's known brightness (K s = 5.09; K s -M' ~ 
0) to calibrate (3 Pic b's brightness. We estimate AM' « 
8.02 ± 0.50 (Mm' « 9.99), where we consider the uncer- 
tainties in our PSF fitting scaling, the dispersion in indi- 
vidual planet magnitude estimates drawn from separate 
reductions, and the intrinsic signal-to-noise of our detec- 
tion. We determine an L' contrast of AL' — 7.71 ± 0.06. 
Combining the L' measurement with the K s band and 
[4.051 data from lBonnefov etaTI (|2011ft and lQuanz et al.l 
( 2010), we have three good quality photometric points to 
investigate the family of possible solutions for (3 Pic b's 
atmospheric properties. 

Figure [4] compares the (3 Pic b photometry to best- 
fit spectra for models with log(g) = 3.5/4-4.5 and T e ff 
= 1000-1800 K for a range of cloud prescriptions: the 
Mode l A and AE thick clou d prescriptions respectively 
from lCurrie et al.l (|2011ft and lMadhusudhan et all ([20111 ) 
that best fit the HR 8799 planet SEDs, the Model E 
cloud deck prescript ion appropriate for brown dwarfs 
(jBurrows et al.l 12006] ) . and a cloudless atmosphere. The 
xl values for these models are, respectively, xl — 24.8, 
12.3, 20.4, and 43.2 for Models A, AE, E, and the cloud- 
less case. The AE thick cloud model provides the best 
fit. Thick cloud models also produce redder L'-M' col- 
ors at high temperatures, similar to that estimated here 
(« -0.22 ± 0.50), though our lack of a reliable M' pho- 
tometric calibration precludes strong conclusions. Good 
photometry is available in only three filters, so we cannot 
yet say that {3 Pic b has thick clouds like the HR 87 99 
planets ([Currie et al.ll2011l ; IMadhusudhan et al.ll201lD . 

The range of gravities and effective temperatures are 
log(g) = 3.5-4.5 and T eff = 1400-1800 K. The implied 
masses for these models range between 4.1 Mj and 19.2 
Mj and ages range between 1 and 27 M yr, broadly con- 
sisten t with the planet mas s (9 ± 3 Mj; ILagrange et al.l 
12010ft . stellar age (12 Myr; iZuckerman et all 12001ft. and 
likely formation timescale (< 3-5 Myr; ICurrie et al.l 
2009). Planet fluxes in the near-IR (1-1.65 /zm) and 
the 3-3.5 /im methane absorption trough are highly 
sensitive to cloud struct ure (e.g. ICurrie et al.l 120111 : 
IMadhusudhan et al.l l20TTft . Thus, J or H broadband 
data and/or narrowband 3-3.5 [im data will be critical 
in breaking model fitting degeneracies. 

4. DISCUSSION 

We present a new detection of /3 Pic b in M' band 
at a separation of ~ 0.21" (a pro j ec t d = 4.05 AU) from 
archival VLT/NaCo data taken in November 2008 and 
a high SNR rereduc tion of L' data first reported by 
ILagrange et al.l (|2010t ). using these data to constrain the 
planet's orbit and atmospheric properties. For orbits 
whose fit to the data yield x 2 < 2.23, we find that the /3 
Pic planet has a semimajor axis of a p = 10.99 AU^'giAU 
and a moderate/low eccentricity (e < 0.31). Admitting 
orbital solutions with x 2 < 8, the parameter ranges are 



9.50 AU 



+3.93AU 
-1.7AU 



and e < 0.23. In both cases, values 



for the planet's inclination (i ~ 88.06-89.69 degrees) and 
longitude of ascending node (SI ~ 30.56-32.12 degrees) 
are tightly constrained and imply that the planet's or- 



bit is almost perfectly aligned with the outer debris disk, 
but not the inclined inner disk (Jl ~ 35-36 degrees) . We 
cannot extract reliable photometry from our M' band 
data; new data at 1-1.65 /xto and ~ 3-3.5 fj,m is needed 
to constrain j3 Pic b's atmosphere. 

Numerous studies of the (3 Pic debris disk(s) have 
identified the star as harb oring a young planetary sys- 
tem fe.g.lSmith and Terrildll984tlKalas and Jewittlll995l: 
iMouillet et al.lll997t iWeinberger et al.ll2003D . More re- 
cently, the presence of a warp in the disk at ~ 80 AU 
- due to the combined effects of the main disk with 
PA ~ 30.8 and a second, inclined disk offset by 5 de- 
grees - was ide ntified as a clear signpost of a perturb- 
ing planet ("e.g. IMouillet et al.l U997T: iHeap et al.l 120001: 
iGolimowski et al.ll2006ft , motiva ting high-contra s t ima; 
ing studies to image the planet. ILagrange et al.1 (2009 
identified (3 Pic b as a likely source of the inclined disk 
and use d the disk morphology to derive mass estimates 
(see also ILagrange et al"1l2010ft . 

Our results suggest that /3 Pic b is probably not aligned 
with the inner disk/warp but rather the main disk, as 
the allowed ran ge in fl is offset f rom t he main disk 
as measured bv IKalas and JewittI (|1995ft by no more 
than ~ 1 degree. Furthermore, the pl anet may be 
misaligned with the submm disk emission (|Wilner et al.l 
12011ft. which is sensitive to dyn amical sculpting by plan- 
ets (|Kuchner and StarU 12010ft . However, models ac- 
counting for the inclined inner disk presume that the 
plan et's orbit is also incl i ned relative to the ma in disk 
(e.g. IMouillet et all Il997t lAugereau etall l2001h . New 
beta Pic b astrometry, a more precise astrometric cali- 
bration of existing beta Pic NaCo data by determining 
and correcting for image distortion, and a detailed rel- 
ative calibration between NaCo data and data reveal- 
ing the disk will further clarify how beta Pic b's or- 
bital plane compares to that for the main disk and the 
inner disk/warp. Furthermore, our new orbital con- 
straints for /3 Pic b strongly motivate new studies of 
the dynamical sculpting of (3 Pic's debris disk by plan- 
ets. If (3 Pic b or non-plan et related mechanisms (e.g. 
IArmitag~ and Pringlc 1997) fail to explain the inclined 
debris disk/warp, the existence of additional planets in 
the system may be required. 

Our M' band detection demonstrates that it is possible 
to directly image planets at separations approaching the 
telescope diffraction limit without sparse aperture mask- 
ing interferometry (SAM: [Ireland and Krausl l2008). The 
high Strehl ratio, large amount of field rotation, large 
mid-IR planet brightness, and LOCI processing pipeline 
are the keys to closing this gap. While SAM can de- 
tect planets interior to the telescope diffraction limit, it 
is overall less sensitive. However, the techniques can be 
complementary, yielding detections or robust limits on 
infant gas giant planets around the youngest stars on ~ 
5-100 AU scales. 

Upcoming facilities like GPI, SPHERE, SCExAO, 
and Project 1640 achieve higher contrast at small in- 
ner working angles in the near-IR p rimarily through 
more sophisticated wavefront control (Macintosh et al.l 
20081: iBeuzit et al.l 120081: iMartinache and Guvonl 120091 : 
Hinklev et al.l 120 111 ). Our results , coupled with previous 
L' ba nd detections of /3 Pic from ILagrange et all ([2009a, 



12010ft and the high signal-to-n oise L' band detection of 
HR 8799e (jMarois et al.ll2010D suggest that the mid-IR 
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TABLE 1 
13 Pic Data Used in this Paper 



Astrometry 

Date 



11-10-2003 

11- 11-2008 

12- 29-2009 
03-20-2010 
Photometry 
Date 



03- 20-2010 
12-29-2009 

04- 03-2010 
11-11-2008 



Filter 

~I7 

M' 

V 

K s 

Filter 
V 

[4.05] 
M' 



Separation ("), Position Angle (°) Reference 



0.411 ± 0.008, 31.7 ± 1.3 
0.210 ± 0.027, 211.49 ± 1.9 
0.326 ± 0.013, 210.64 ± 1.2 
0.345 ± 0.012, 209.8 ± 0.8 

Absolute Magnitude 
11.20 ± 0.12 
9.73 ± 0.06 
9.77 ± 0.23 
« 9.99 ± 0.50 



Lagrange ct al. (2009a) 
this work 
this work 

Bon nefov et al. (2011) 



Reference 
Bonnefov et al. 



(2011) 

this work 

Quanz et al. (2010) 
this work 



Note. — Our astrometry and photometry are drawn from the three separate 
reductions shown in Figure 1. The L' measurement assumes L^ pic = 3.45. The M' 
photometry lacks reliable photometric calibration and thus is not useful for atmospheric 
modeling. 



may also be fertile ground for new exoplanet detections at 
small separations for very young systems. Young, nearby 
1.5-2 Mq stars like /? Pic are partic ularly promising tar- 
gets for direct imaging surveys (e.g. iCrepp and Johnson! 
I20T1 and many have resolved debris disks (e .g. HD 
181327. ISchneider et al. ll2006l:IChen et~aT~l[2008h . Imag- 
ing massive planets in such systems can yield additional 
studies of planet-disk interactions, such as those moti- 
vated by this work. 
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Fig. 1. — (Top-left) Reduced VLT/NaCo L' band image from Dec. 29, 2009 data showing a 21-ct detection of /3 Pic b at a separation 
of ~ 0.32". (Top-right) Reduced image showing the detection of ft Pic b in M 1 band (same image size). The LOCI parameters used to 
construct the reduced image include 5 = 0.25, g = 0.4, and Na = 250xFWHM. (Bottom panels) Reduced images using less aggressive 
LOCI settings - S = 0.5, g=0.5, and N A = 300xFWHM (left panel) and S = 0.5, g=l, and N A = 300xFWHM (right panel) - yielding 
detections of SNR ~ 5 and 4.5, respectively. In general, we detect the j3 Pic planet at a 4-6 a level using a range of LOCI parameters: 8 
= 0.25-0.5, g = 0.3-1, N A = 200-300. 
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Fig. 2. — Orbital analysis results following the Monte Carlo method outlined in Thalmann ct al. (2009), showing parameters for solutions 
fulfilling x 2 < 2.23. The left panel displays the family of solutions in semimajor axis, eccentricity, and inclination space. The right panels 
show histogram distributions of these three parameters and the longitude of the ascending node. The vertical dashed line identifies the 
weighted median value for each parameter; the vertical dotted lines define the 68% confidence interval for each parameter. 
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Fig. 3. — Same as previous figure except for orbital solutions fulfilling \ < 




Fig. 4. — Figures comparing the /3 Pic b photometry to spectra for a range of cloud prescriptions : the Model A 'thick cloud' limit used in 
ICurrie ct al. (201l[) (upper-lef t), the Model AE 't hick cloud' prescription from Madhusudhan ct al. (201f ) (uppe r-right) , the M odel E brown 



dwarf cloud prescription from[Burrows ct al. (2006) (lower-left), and a cloud-free atmosphere also from Burrows ct al. (2006) (lower-right). 
The blue dots show the flux predictions for each bandpass. All models assume a modal particle size of 60 fim. 



